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Abstract 

Atom interferometers can be used to study phenomena leading to ir- 
reversibility and dissipation, induced by the dynamics of fundamental 
objects (strings and branes) at a large mass scale. Using an effective, 
but physically consistent description in terms of a master equation of 
Lindblad form, the modifications of the interferometric pattern induced 
by the new phenomena are analyzed in detail. We find that present 
experimental devices can in principle provide stringent bounds on the 
new effects. 
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1. INTRODUCTION 



The evolution in time of a system S immersed in a large environment £ can be obtained 
from the dynamics of the total system S + £ by eliminating (i.e. integrating over) the 
degrees of freedom of £. Since S + £ is closed, the total dynamics is unitary; this is no 
longer true for the evolution of the subsystem S alone, which in general turns out to be very 
involved, developing non-linearities and memory effects. However, when the interaction 
between subsystem and environment is weak and there are no initial correlations between 
S and £, the time evolution of S can still be realized through linear maps on the states of 
iS, satisfying basic physical requirements, like forward in time composition law (semigroup 
property), entropy increase (irreversibility) and complete positivity (that guarantees the 
physical consistency of the evolution in all situations). These one-parameter (=time) 
family of maps form a so called quantum dynamical semigroup [1-4] and are generated by 
a master equation of Lindblad form [5] . 

This description of the time-evolution of open systems is very general; it was originally 
developed in the framework of quantum optics [6-8] , and subsequently used to model very 
different physical situations, from the study of various statistical systems [1-3], to the 
analysis of the interaction of a microsystem with a macroscopic measuring apparatus [9- 
11], to the description of the emergence of the classical world [12, 13] and of the so called 
dynamical reduction [14]. 

Master equation of Lindblad form can also be used to describe phenomena leading 
to irreversibility and dissipation at low energy generated by the dynamics of fundamental 
objects at a large scale, typically the Planck mass [15]. Indeed, the dynamics of extended 
objects, strings and branes, gives rise at low energies to a weakly coupled heat bath, and 
as a consequence to decoherence phenomena [16]. From a more phenomenological point of 
view, similar effects have also been described in the framework of quantum gravity [17]: 
due to the quantum fluctuations of the gravitational field and the possible generation of 
virtual black holes, spacetime should become "foamy" at scales comparable to the Planck 
length, leading to loss of quantum coherence [18-23]. Furthermore, dissipation and deco- 
herence are also the natural outcome of the general dynamics in theories with large extra 
dimensions [24]: indeed, the possible energy leakage from the boundary of spacetime (our 
four-dimensional brane universe) into the bulk due to gravity effects would inevitably inject 
noise into the boundary, thus inducing irreversibility and dissipation at low energy. 

Our present knowledge of string theory does not allow precise estimates of the magni- 
tude of these non-standard effects. Nevertheless, dimensional arguments suggest that they 
must be very small, being suppressed by at least one inverse power of a large fundamental 
mass. Despite of this, they can be in the reach of various interferometric devices. Indeed, 
detailed investigations involving different elementary particle systems (neutral mesons [25, 
26], neutrons [27], photons [28], and neutrinos [29]) have shown that present and future 
experiments might soon reach the sensitivity required to detect the new, non-standard 
phenomena. 

Another physical situation in which phenomena leading to irreversibility and dissipa- 
tion can be studied is provided by atom interferometers, where a beam of nearly monoen- 
ergetic atoms is coherently split into two components that are recombined at the exit of 
the apparatus [30-33]. The interferometric pattern observed at the end of the devise is 
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influenced by the action of external phenomena, produced e.g. by an external electric or 
magnetic field, or by earth gravity. The sophistication of present interferometric apparata 
is so high that the theoretically predicted changements in the interferometric figure for 
some of these phenomena have been confirmed with high accuracy [30-32]. 

Irreversibility and dissipation also affect the propagation of the atoms in the inter- 
ferometer; this leads to a deformation of the corresponding interferometric pattern at the 
exit of the apparatus. It turns out that these modifications are very distinctive of the 
dissipative phenomena, and can not be mimicked by other physical effects, as the ones 
mentioned before. 

In the following, we shall analyze in detail these modifications under the hypothesis 
that the generalized dynamics of the atoms inside the interferometer be generated by a 
master equation of Lindblad form. For sake of definiteness, we shall limit our discussion 
to three-grating atom interferometers in the Bragg regime, where the split and the re- 
combination of the incident beam is realized by material gratings or laser standing waves. 
The approach has points in common with the one adopted in [27] to study similar effects 
in neutron interferometry [34]. The perturbative treatment adopted there is however in- 
applicable in the present case: this fact, together with the operational differences in the 
actual functioning of an atom interferometer require a completely different and indepen- 
dent analysis. As discussed in the final section, the outcome of our investigations is that 
atom interferometry experiments could provide the most accurate estimate of the non- 
standard, dissipative effects that can be induced by a fundamental dynamics at Planck's 
scale. 



2. MASTER EQUATION 

The evolution of the atoms inside the interferometer can be analyzed using an abstract, 
two-dimensional Hilbert space. The states corresponding to the two split beams in the 
apparatus can be taken to be the basis states in this space. More in general, the quantum 
state of an atom travelling inside the interferometer will be a statistical mixture of the 
basis states, and therefore described by a density matrix p, i.e. by a hermitian, positive 
defined operator with unit traced With respect to the chosen basis, one can then write: 

P=( P n 1 P n 3 ) , P4=P3, Pl+P2 = l, (2.1) 

V Pi J 

where * signifies complex conjugation. 

As explained in the introductory remarks, the starting point of our analysis is the 
assumption that the dynamics of the atoms inside the interferometer be generated by a 
master equation of Lindblad form 1 '- [1-5] : 

^ t p(t) = -i[H lP (t)}+D[p(t)] . (2.2) 

t For earlier works on the use of the formalism of density matrices in atom interferom- 
etry, see [35-37, 33] and references therein. 

* An equation of this type has also been used to study decoherence effects in posi- 
tion space induced by the scattering of photons on the atoms inside the interferometer 
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The first term in the r.h.s. represents the standard hamiltonian contribution. In the chosen 
basis, the effective hamiltonian can be written as: 



H = 



E + u) 








E 



(2.3) 



where E is the energy of the atoms in the incident beam. On the other hand the splitting 
in energy 2ui among the two internal beams is usually induced by the action of laboratory 
controlled effects, typically the presence of external fields. For open system though, even 
in absence of external fields, the quantity u> is in general nonvanishing. Indeed, one can 
show that the weak interaction of the system with the external environment induces in 
general a hamiltonian contribution, giving rise to the "Lamb shift" term u in (2.3) [1-3, 
16, 29]. 

Nevertheless, it is the additional piece D[p] in (2.2) that describes true mixing en- 
hancing phenomena: in absence of it, the evolution of p would be unitary and reversible. 
It can be represented by a trace-preserving linear map acting on the three independent 
components of the density matrix in (2.1). Decomposing for convenience ps in its real and 
imaginary parts, 

Ps = P 1 ~ ip 2 , (2.4a) 

and introducing the combination 

Pl -p 2 = 2p 3 , (2.46) 

one can then write D[p] as a 3 x 3 real, symmetric matrix V, acting on the real vector \p) 
of components (p 1 , p 2 , p s ): 

a b c 

V = -2 b a (3 . (2.5) 
_c 7_ 

The six parameters a, 6, c, a, (3 and 7, with a, a and 7 non negative, are not all independent; 
physical consistency of the full time evolution (i.e. the request of complete positivity, see 
[38, 39] for details) further imposes the following inequalities: 

RS > b 2 , 
RT> c 2 , 



2R = a + -/-a>0, 
2S = a + -f-a>0, 
2T = a + a- 7 >0, ST > (3 2 

RST >2bcf3 + Rf3 2 + Sc 2 + Tb 2 . 



If one includes also the hamiltonian contribution and further recalls that Tr[p(f)] = 1, 
the evolution equation (2.2) can be rewritten as a diffusion equation for the 3-vector \p(t)): 



d_ 
dt 



\p(t)) = -2H\p(t)) , H 



a b + u c 

b — uj a (3 
c Pi 



(2.7) 



[37]. Instead, the evolution equation (2.2) is written in "polarization" space, and, as ex- 
plained below, represents the most general master equation compatible with basic physical 
requirements. 
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Its solution involves the exponentiation of the matrix 7i: 



\p(t))=M(t)\p(0)) , M{t) = e~ 2H \ (2.8) 

where |p(0)) represents the initial state of the atoms entering the interferometer. It coin- 
cides with one of the following density matrices 

corresponding to the possible choices of orientation of the incident atomic beam with 
respect to the first diffracting grating. Both choices lead to the same final results; for 
definiteness, in the following we shall work with p^\ so that |p(0)> = (1/2,0,0). 



3. OBSERVABLES 

In the language of density matrices, physical observables are represented by suitable 
hermitian operators, whose mean values can be obtained by taking their trace with p(t). 
In particular, the intensity pattern observed at the end of the interferometer is given by 
the mean value of the following projector operators [18, 27]: 

that correspond to the two possible exit beams in which an atom can be found while exiting 
the apparatus. In the standard situation, it is the phase 9 that gives the modulation of the 
interferometric pattern. This is usually obtained by moving the transverse position of one 
of the gratings (or laser standing waves) responsible for the diffraction of the atom beam. 
Indeed, in an idealized situation one finds [30-32]: 

6 = k(x 1 -2x 2 +x 3 ) , (3.2) 

where Xi, i = 1, 2, 3, represents the position, transverse with respect to the incident beam, 
of the z-th grating, while k is the wave vector of the diffracting lattice. 

The intensity I± of the interference figure detected at the two possible exits is then 
given by: 

I±(t) = (0±) = Tr[0 ± pit)] = \ + O x ± p\t) + 0\ p\t) + Ol p 3 (t) , (3.3) 

where definitions similar to the ones in (2.4) have been introduced also for the entries of 
the two matrices 0±. Using (3.1), one finds: 



J ± (t) = I|l±2 cos 9 p 1 (t) + sin 9 p 2 (t) j. (3.4) 
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Since an atom exiting the interferometer can only be found in one of the two exit beams, 
particle conservation requires I+(t) +Z_(t) = 1, which is clearly satisfied by (3.4). 

The intensity curves in (3.4) can be compared with the experiment, provided explicit 
expressions for the entries of the matrix A4(t) in (2.8) are given. Formally, this can be 
obtained by studying the eigenvalue problem for the 3x3 matrix 7i in (2.7): 

H \v k) = Afc \vk) , k= 1,2,3. (3.5) 

The three eigenvalues Ai, A 2 , A 3 satisfy a cubic equation, 

X s +rX 2 + s\ + w= 0, (3.6) 

with real coefficients: r = — (a + a + 7), s = aa + + 07 — b 2 — c 2 — (3 2 + u 2 , w = —detH. 
It then follows that A& are either all real, or one is real and the remaining two are complex 
conjugate; further, in both situations, one can show that in presence of dissipation the 
three eigenvalues have always positive real parts [40]. 

Using the fact that the matrix H itself obeys the equation (3.6), one finds: 



1 f P -2Afct r, x 

I±(<) = 2 { 1 ± E 3AfT2^T7 ( A ' - {a + ^ + <" - 13 ) C ° S 6 

k=l K 

+ ((a;-6)(A fc -7)- / 9c) sin^ | . 



(3.7) 



From this general expression, one sees that in presence of complex eigenvalues, a further 
harmonic modulation in time of the interference figures occurs, while exponential damping 
terms always prevail for long enough times. Further, note that in absence of dissipation, 
a = b = c = a = f3 = "f = uj = 0, and thus = 0, the expressions of T± in (3.7) reduce 
to their standard, time-independent ones [30-32]: 

J± = ^{l±cos#} ; (3.8) 

any deviation from this formula as described by (3.7) clearly signals the presence of dissi- 
pative phenomena in atom interferometry. 

Although explicit expressions for the eigenvalues A^ can always be found via Cardano's 
formula [41], the form (3.7) of the intensities I± is rather involved, and of limited use in 
practice. Having in mind possible comparison with experimental data, the study of suitable 
approximations of (3.7) might result appropriate. 

In this respect, a useful working assumption is to take 7 = 0;^ in this case, the 
inequalities (2.6) further impose b = c = (3 = and a = a. In this simplified situation, 



t There are essentially two known ways of implementing the condition of weak inter- 
action between subsystem and environment [1-3]: the singular coupling limit (in which 
the time-correlations in the environment are assumed to be much smaller than the typical 
time scale of the subsystem) and the weak coupling limit (in which it is the subsystem 
characteristic time scale that becomes large). One can check that the second situation 
leads precisely to the condition 7=0 [29] . 
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the formula in (3.7) reduces to: 



J ± (t) = 1 1 1 ± e ~ 2at cos (9 - 2wt) } . (3.9) 

This is surely the most simple expression that the intensity probabilities T± (t) take in pres- 
ence of dissipative effects. It differs from the standard expression in (3.8) by the presence 
of an exponential damping factor and of an additional harmonic phase that accumulates 
in time. 

A different approximation of the general formula (3.7) can be obtained when the 
parameters a, 6, c, a, (3 and 7 are small with respect of u. This could happen when the 
interferometer is immersed in a strong external field, so that its contribution to the energy 
shift u) due to its interaction with the atom beams largely overrides the one coming from 
the effects of a weakly coupled environment. * In this case, the additional piece D[p] in 
the evolution equation (2.2) can be treated as a perturbation. Using the solution of this 
equation expanded up to second order in the small parameters, from (3.4) one obtains: 

Z±{t) = X - /lie" ( a+a >* fcos2Qt + ^^sin2fit- ^ sin 2 fit") cos6> 

/&-« c0 \ (3 ' 10) 

+ ( — — sin2fit- ^-cos2fit ) sinfl 
\ il \V J 

where O = [u 2 — b 2 — c 2 — f3 2 — (a — a) 2 /4] 1 / 2 . In writing (3.10) we have reconstructed 
the exponential factor by putting together the terms linear and quadratic in t; a similar 
treatment has allowed writing all harmonic pieces in terms of the frequency O. It is worth 
noting that for c = (3 = the formula (3.10) gives the exact expression for the intensities 
T±: no approximation is involved. This is a consequence of the fact that for c = (3 = the 
matrix 7i in (2.7) becomes block diagonal, so that explicit, manageable expressions for its 
exponential Ai(t), and therefore for I±, can be given. From this point of view, the validity 
of (3.10) goes beyond the second order approximation in which it has been derived: it can 
be considered as the expansion of the full expression (3.7) for Z± up to second order in c 
and f3. 



4. INTERFERENCE PATTERN 

The behaviour of the general expression (3.7) for X± and of its special cases (3.9) and 
(3.10) crucially depend on the dissipative parameters: at least in principle, they can be used 
to obtain informations on their values from fits with the experimental data. The magnitude 
of o, 6, c, a, (3, 7 and u are nevertheless expected to be very small. In fact, for subsystems 

* In this respect, it should be noted that even in absence of any external field, a hierarchy 
between u and the other dissipative parameters a, 6, c, a, (3 and 7 could be nevertheless 
generated by the interaction with the environment. For details, see [16] and the Appendix 
in [29]. 
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in interactions with large environments, the effects leading to dissipation and decoherence 
can be roughly estimated to be proportional to the typical energy of the system, while 
suppressed by inverse powers of the characteristic energy scale of the environment [1-3]. 
In the case of non-standard phenomena induced by the dynamics of fundamental objects 
(strings, branes) at Planck's mass Mp, an upper bound on the magnitude of the dissipative 
parameters can be roughly evaluated to be of order M\/Mp, where Ma is the mass of the 
atoms used in the interferometer [16, 29]; in typical real situations, this ratio takes values 
between 10 -18 and 10 -15 GeV (or equivalently, between 10 3 and 10 6 KHz). 

A further difficulty in comparing the theoretically predicted interference figures with 
the experimental data arises from the fact that the previously derived expressions for T± 
hold in the case of an idealized interferometer, with perfectly monoenergetic atomic inci- 
dent beams. In practice, the values of the atom momenta spread over a finite distribution. 
This fact, together with the inevitable imperfections in the construction of the actual 
interferometric apparatus, produce attenuation in the intensity of the signal. 

One can take into account these spurious effects by modifying the previously derived 
intensity spectra with the introduction of suitable unknown parameters. To keep the 
discussion as simple as possible, we shall concentrate on the expression (3.9) for I±; similar 
arguments apply to the other formulas. By denoting with J\f± the atom countings at the 
two exit beams of the interferometer, one generalizes the spectra in (3.9) as: 

A/±(t) = A/± 0) |l ± C± e~ 2at cos (9 - 2ut) j . (4.1) 

The constants C± are the fringe contrast and parametrize the intensity attenuation, while 
A/l 0) are suitable normalization factors;* note that particle conservation now requires: 

A/"|°' ) = A/"i°' ) C-. Clearly, the higher the fringe contrast, the more accurate the deter- 
mination of the dissipative parameters a and u> from the experiment will be. 

In order to fit actual experimental data with the expression (4.1), further elaborations 
are however needed. As mentioned before, the intensity spectra are reconstructed by 
counting the atoms at one of the exit beams as a function of the transverse position x of 
the final grating (or standing laser wave), with respect to a reference, initial situation. This 
means that the geometry of the two paths followed by the atoms inside the interferometer 
slightly changes as x varies;* as a consequence, also the total evolution time t, the time 
spent by the atoms inside the interferometer, changes with x. Since the two path inside 
the apparatus are very close to each other (their actual separation is at most 100 //m), for 
small x one finds: 

#1 

t = t + — x , (4.2) 

v 

* In absence of dissipative effects, a theoretical estimate of C± has been obtained using 
atom optics [42, 43]. 

* The situation is completely different in a neutron interferometer [34]: made of a 
silicon crystal, its geometry can not be varied. In this case, the interferometric spectra are 
obtained through a thin slab of material inserted transversally to the two beams inside 
the interferometer; a slight rotation of it produces a phase difference u between the two 
"optical" paths. 
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where v is the average velocity of the atoms in the incident beam, while #1 is the first 
order Bragg diffraction angle (typically of order 10 -4 rad). On the other hand to is the 
fixed time of flight of the atoms inside the interferometer when it is in its initial, reference 
status; it can be determined with high accuracy from the geometric specifications of the 
actual apparatus and can be modified only by changing the energy of the primary atom 
beam, or by modifying the longitudinal dimension of the interferometer. 

The outcome of this discussion is that it should be possible to estimate the values 
of a and 10 from the behaviour in x of the expression in (4.1), and therefore from a fit 
with experimental data. Unfortunately, in the experimental set-ups so far constructed the 
dependence of t on x can hardly be seen: one finds that while to is at most 10~ 3 sec, the 
quantity "&\x/v results at least ten orders of magnitude smaller, even for maximal values 
of the displacement x (a few hundreds nm). 

Therefore, as a good approximation one can safely take t ~ to, and rewrite the inter- 
ference pattern (4.1) as 

Af±(x) = A/l 0) {l ± [p± cos(6> + kx) + Q± sin(0 o + } , (4.3) 

where 

V± = C± e~ 2at ° cos 2ut , Q± = C± e~ 2at ° sin 2ut , (4.4) 

while #o is a fixed phase that is characteristic of each interferometer. A fit of (4.3) with 

experimental data will allow to determine the parameters a/1 o) , V±. l Q± and 6q, and 
therefore to obtain informations on the dissipative parameters. 



5. DISCUSSION 

We shall now briefly report on the results of a x 2 fit of the formula (4.3) with recently 
published data from two experiments [44, 45]. The analysis that follows is of limited 
quantitative meaning: direct access to the data and a careful study of systematic errors 
are needed in order to obtain precise determination of the dissipative effects; nevertheless, 
it will provide a rough estimate about the sensitivity of present atom interferometers to 
the parameters a and to. 

The atom interferometers used in the two experiments are particularly sensitive de- 
vices, reaching a very high fringe contrast; this is obtained by using neon, respectively 
lithium, atoms as "matter waves" and laser standing light as diffracting device. Both in 
[44] and [45], only data from one of the two exit beams are reported, the ones corresponding 
to the lower sign in (4.3). Since as mentioned before to is known with high precision, from 
the ratio V-j Q- one can immediately obtain an estimate for the parameter ui. Using the 
data from the first experiment, from our fit we find uj = (0.7 ± 0.2) x 10~ 21 GeV, where 
the quoted error is only statistical. 

On the other hand, the determination of a is subordinated to the estimate of the 
fringe contrast C_. This would not be necessary if the parameters V- and Q_ can be 
measured for two different values of the flight time to- As mentioned before, this can be 
obtained either by changing the average velocity of the incoming atoms, or by varying the 
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dimensions of the interferometer. In lacking of this extra information, we shall obtain an 
estimate for C_ using directly the data. 

In absence of dissipative effects, a = uj = 0, the constant C_ can be obtained from 
the maximum, 7V (max) , and the minimum, , atom counts in the experimental inter- 

ference figure: C_ = (A/" (max) - A/"( min ))/(A/" (max) +W (min) ). Although this formula is only 
approximately valid for nonvanishing a and uj, in practice the systematic error that one 
makes in adopting it can be estimated to be at the end much smaller than the pure ex- 
perimental uncertainty. Using the rough experimental data, one then deduces: C- ~ 62%. 
With this value, one finally gets: a = (0.1 ± 0.1) x 10 -22 GeV, which is compatible with 
zero. 

The accuracy in the determination of a and uj improves using the data from the most 
recent experiment, thanks to the higher fringe contrast (of about 74%) and the increase in 
the number of experimental points. In fact, the same procedure adopted before now gives 
the following estimates: a = (0.3 ±0.1) x 10" 23 GeV and uj = (0.20 ±0.01) x 10" 21 GeV. 1 
Note that these values are perfectly compatible with the ones previously determined in 
the case of the neon interferometer. As explained before, the values of the dissipative 
parameters should be proportional to the square of the mass of the atoms in the incident 
beams. Therefore, the values of a and uj determined with the data from the lithium beam 
should result smaller than those obtained from the first experiment. 

In conclusion, the results of our discussion show that atom interferometers are poten- 
tially very sensitive to the presence of phenomena leading to dissipation and decoherence. 
Although the performed error analysis has been limited to statistical uncertainties, the 
derived estimates seem to indicate nonvanishing values for a and u, of magnitude compat- 
ible with an origin from a fundamental dynamics at a very large mass scale. As already 
remarked, direct access to the rough experimental data and more complete x 2 fits are 
needed in order to claim the presence of dissipative effects. We nevertheless hope that our 
preliminary analysis will stimulate further, more accurate investigations. 



t We remark that as before the quoted errors are purely statistical; a thorough analysis 
of the full experimental data, that takes into account also the systematics, would likely 
worsen the estimated errors, in particular that on uj. 
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